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Abstract 


Using the nonequilibrium Green function, we show that microwave irradiation 
can reverse the supercurrent flowing through a superconductor / quantum- 
dot / superconductor structure. In contrast with the conventional sideband 


effect in normal-metal / quantum-dot / normal-metal junctions, the photon- 


assisted structures appear near Eo = $hw (n = +1,+2---), where Ep is the 


resonant energy level of the quantum dot and w is the frequency of microwave 
field. Each photon-assisted structure is composed of a negative and a positive 
peak, with an abrupt jump from the negative peak to the positive peak around 
Eo = .سي‎ The microwave-induced m-junction transition is interpreted in the 
picture of photon-assisted Andreev bound states, which are formed due to 
multiple photon-assisted Andreev reflection between the two superconductors. 


Moreover, the main resonance located at Eg = 0 can also be reversed with 


proper microwave strength and frequency. 


PACS numbers: 74.50.+r, 73.63.Kv, 72.23.Ad. 

When two superconductors are weakly linked, dc current can flow even without bias 
voltage. The driven force of the supercurrent is the phase gradient in the macroscopic wave 
function of Cooper pair condensate. The supercurrent and the phase difference across the 
junction has the relation J = /,sing, with I. > 0 being the critical supercurrent. If the 
link area of the junction is controlled by some external conditions, J, may be enhanced , 
suppressed, or even reversed. The reverse sign of 7, is referred as to m-junction transition, 
since the minus sign can be absorbed into ¢ as an internal phase shift of 7. The ground 
state of m-junction is ع = و‎ rather than o = 0 as in usual 0-junction. This is easily seen 
from the relation / = mm and therefore F — —#I, cos @+ Fo, where F is the free energy of 
the junction. The studies of m-junction have not only academic interests, but also potential 
applications, e.g., realization of qubit in a superconducting loop with 0- and m- junctions 
uj. 

In addition to the intrinsic 7 phase shift in the order parameter of high Te superconduc- 
tors due to the d-wave symmetry [J], several mechanisms were also proposed and investigated 
for carrying out m-junction with the conventional BCS superconductors. (1) Coupling two 
superconductors through a ferromagnetic layer: 'The idea can be traced back to the original 
works of Fulde and Ferrel [B] and Larkin and Ovchinnikov [fH], who independently predicted 
that superconducting order parameter can be modulated by an exchange field, and contains 
nodes where the internal phase shift is 7. With the growing interests in the hybrid system 
of superconducting / ferromagnetic (S/F) materials, m-junction behavior in SFS sandwiches 
and SF superlattices were intensively studied [E48]. Very recently, experiments on Joseph- 
son junction inserted with a weakly ferromagnetic interlayer did observe the m-junction 
transition [D]. (2) Coupling two superconductors by an Anderson impurity or a quantum 
dot: The early papers of Glazman and Matveev and Spivak and Kivelson showed 
that when the impurity is singly occupied due to Coulomb repulsion, the sign of Josephson 
current is opposite to that without the repulsion.. In the last decade, a number of theoret- 
ical papers were devoted to this issue [12] [I5]. Experimentally, the technique of fabricating 
Josephson junctions containing nanoparticles was available [[6], yet no relevant results on 
T-junction transition was reported. (3) Introducing a nonequilibrium distribution in the cen- 
tral mesoscopic region: In a mesoscopic superconductor / normal-metal / superconductor 


(SNS) structure, the quasiparticle distribution can be driven far from the equilibrium by a 


control voltage across the N region. When the control voltage exceeds a certain value, the 
nonequilibrium distribution has so much weight on the negative part of the current carrying 
density of states that the supercurrent reverses its sign. After a few years of the prediction 
[17-20], the reverse of the supercurrent was successfully observed in a controllable Josephson 
junction BIJ. 

In this paper, we propose a new mechanism for the realization of m-junction, schemat- 
ically shown in the inset of Fig.1. Consider a Josephson junction consist of two BCS su- 
perconductors (S) coupled through a quantum dot (QD). Applying a microwave (MW) on 
the QD, the resonant energy level of QD will shift adiabatically with the time dependent 
external field. We shall show below that 7-junction transition occurs in the photon-assisted 
Josephson current, and the main resonance can also be reversed by a proper choice of MW 
strength and frequency. The MW-induced z-junction transition is related to the formation 
of photon-assisted Andreev bound states (PAABS), which is a generalization of the usual 
Andreev bound states (ABS). 

We model the S-QD-S system by the following Hamiltonian, 


H = H; + Hp + Hr + Hr , (1) 


where Hg = DE OT Uke tX | جره بره تسوه‎ + H.c. with 8 = L / R is the 
standard BCS Hamiltonian for the left / right superconducting electrode; Hp = (Eo + 
W coswt) .رر‎ cle, is for the QD under MW irradiation, in which the intradot interaction 
is ignored for simplicity 22], and Hr = 3) Bas (CBG io + H.c.) is the tunnel Hamiltonian , 
connecting the three parts together. 

The time dependent current flowing through S-QD-S can be formulated using the 


nonequilibrium Green function [3], 


Ig(t) = z2 Re in a. | ats IG" (t, 1) UG (t1, t) + G*(t, mast] } (2) 


in which the full Green function of QD and the self energy induced by the coupling with 8 


electrode are defined as 
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and يه‎ in Eq.(2) is the third Pauli matrix, with the diagonal element +1 for the electron 
current (spin?) and -1 for the hole current (spin|). G” and GS obey the corresponding 
Dyson equation and Keldysh equation: 


6" (t3, t2) = g' (t, t2) + f [utei og e e , (6) 
8َ I 1 dtdt G'(t,£)9 (t, C) G"(I ts) . (7) 


The remaining task is to solve these equations properly. It should be pointed out that finite 
perturbation expansion of Eq.(6) is inadequate in the problem, because the formation of 
PAABS involves up to infinite order of tunneling processes. 

Although G” (t1, t2) is no longer the function of tı — t», it still holds that G” (tı + 2 to + 
27) = G" (tı, t9) due to the periodical time dependence of the MW. Hence G” (t1, t2) can be 


Fourier expanded as 


G" (ti, t2) 2 E. e ic - t2) 0) ). (8) 


Define the Fourier transformed Green function 
7ہ‎ (e = G5 (e nw) , (9) 


which has the property that if C(t;,t3) =  [dtA(t,t)B(t,t9) then Crp(e) = 
ر2‎ Amk(€)Bun(€). The Fourier transformed g'and 3^* can be obtained as 


Jema) E RM 0 
ارول‎ - 7 ۱ , _ (10) 
0 itd) cor dnas UU 


Xm) = 30 (e — mu), (11) 


in which J,(o) is the nth Bessel function with the argument a = Z, X*(e) = f(e)[Z*(e) — 
>" (e)] with f(e) being the Fermi function, and X” (ع)‎ = X (e) + Up(e) with 


Mf, = اھ‎ Lo gae” 
8 6 (e+ ig)? — A? =A وفنا‎ 1 


(Im Vz > 0) . (12) 


In the self energy, I'g is the coupling strength between QD and Û electrode, in is the dephas- 
ing rate in the S electrodes which determines the broadening of PAABS. In the numerical 
calculation, the limit 7 — 0* is achieved by extrapolation. 


The Fourier transformed Dyson equation can be expanded as 


G'(c) = g"(e) + g'()2" teg" ) ل‎ + g'(0X' (9g OU (9g (e) +۰... . (13) 


To re-sum up the series, we adopt the resonant approximation: [24 
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which is justified for the weak coupling case of l'a > w. Under this approximation, the 


Dyson equation can be exactly solved as 
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The time dependent current is Fourier transformed as 
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With the relations G<(e) = G"(e)X“(e)G“(e), G^(e) = [G"(e]', and ©") solved in 
Eq.(16), the current can be evaluated. In this work, we concentrate on the time aver- 
aged current I = ny = -1. Notice that I = (Tal — Tai) / (Ur, +I), and de- 
fine © = (UTgX,—T;Eg)/(Tri-4 Tg). Again, we adopt the resonant approximation in 
G<(e) = G'(e)X:*(e)G^(e), and the current formula is simplified as 


7 سے‎ f > h(a) J_y(a) 2Re Tr oz OOF + GS (QE) X (23) 


where G7, (€) = (0) — zr (e) and Ge) = Gy (e) (e)G (c). One can see in the 
formula that / is contributed by various photon-assisted Andreev reflections, in which an 
electron absorbs / photons and a hole emits I’ photons. Eq.(23) is the central result of this 
paper. 

Below we shall numerically study the influence of the MW field on the Josephson current 
flowing through S-QD-S. In the calculation, we set e = h = 1, take kgT' = 0, and choose 
A = 1 as the energy scale. In a qualitative discussion, we regard I. = 1(¢ = 2) as the critical 
supercurrent, and use the sign of I. as the criterion of the 0-junction vs the m-junction. Fig.1 
shows the curves of the critical supercurrent J, vs the resonant energy level Eo, for different 
MW strengths. There are two remarkable features in the plot: (1) With the increase of 
MW strength, the main resonance located at رط‎ = 0 is gradually suppressed, while several 
photon-assisted tunneling structures grow up around Ep = $nhw (n = +1, +2 - - -). (2) Each 
of the structures contains a negative and a positive peak, and J, reverses its sign abruptly 
near Eo = snhw. Meanwhile, the magnitudes of these structures exhibit a non-monotonous 
dependence on the MW strength. Feature (1) is quite different from the sideband effect 
in N-QD-N, where photon-assisted peaks appear near Ey = nhw. In N-QD-S junctions, 
however, photon-assisted peaks also appear near Ey = snhw, which can be explained in 
the picture of photon-assisted Andreev reflection (PAAR) [24]. Basically, feature (1) can 
be understood in the similar way as in N-QD-S, but should be noticed that in S-QD-S the 
photon-assisted electrons and holes undergo up to infinite order of PAAR processes, so that 
PAABS are formed in the QD (schematically shown in the inset of Fig.1). Feature (2) is 
dramatically distinguished from those in either N-QD-N or N-QD-S, indicating that MW 
irradiation can lead to m-junction transition in the photon-assisted current through S-QD-S. 

To understand the feature (2) requires the knowledge of PAABS . For this purpose, 
the time averaged current is rewritten in the form 7 = £f “f(e)j(e). The analysis of 


the current carrying spectrum j(e) provides the information of the supercurrent carried by 


6 


each of the PAABS. Fig.2 illustrates the spectrum of j(¢) for MW strength a = 0.6, with 
Eo = 0, 0.05, 0.10, 0.15, 0.20. For comparison, the spectrum of Ey = 0 and Ep = 0.2 
without MW irradiation are also shown on top of Fig.2. Without MW irradiation, the 
spectrum of S-QD-S has two ABS within the superconducting gap, located near e = +F0, 
carrying supercurrent with opposite signs. The spectrum also has continuous parts outside 
the superconducting gap, negative for w < —A and positive for w > A. Turning on the MW, 
the original two ABS split into two sets of PAABS, due to various photon-assisted processes. 
The PAABS marked with A, locate near e = Eo + nw, carrying supercurrent with positive 
signs; while those marked with B, locate near e = —E + nw, carrying supercurrent with 
negative signs. On the contrast, MW has negligible effect on the continuous spectrum, 
because the contribution from these parts depends strongly on the density of states in the 
S electrodes and hence photon-assisted processes is largely suppressed. Notice that at zero 
temperature J is contributed by the spectrum in the range of e < 0. For Ey = 0 (Fig.2a), 
the states A,, and B,, are energy degenerate, carrying opposite supercurrent weighted by a 
distribution factor J2(o), similar to the sideband in N-QD-N. The contributions from ہے‎ 
and B. ,, A. » and ,وط‎ etc. cancel with each other exactly, and the positive contribution 
from Ap overwhelms the negative contribution from C7, resulting in a positive current. With 
the change of the resonant energy level E, An and رط‎ move toward —A and +A, respectively. 
For Ey = ںا‎ (Fig.2d), the states A, and B,,_; are energy degenerate, and therefore strongly 
hybridized with each other. So the weight of A, and B,,..; are re-distributed, such that they 
carry the same amount of supercurrent with opposite signs. The contributions from Ay 
and B. 4, A-1 and B_», etc. cancel with each other exactly, leaving the negative continuous 
spectrum C7 contribute to the supercurrent. This is the origin of the 7-junction transition in 
the photon-assisted supercurrent. Further investigation shows that the evolution of PAABS 
from (a) to (d) is more complex than stated above. In fact, the hybridization of A, and B, 
already occurs before the energy degeneracy (see Fig.2d and its inset). The hybridizing of 
PAABS obeys the following rule: when two states A; and B; approaches each other, they 
“interact” with each other only if the energy of A; is higher than that of B;. With this rule, 
the abrupt sign change in the current near Ey = hw is readily understood. 

In the above, we have shown that MW irradiation can reverse the sign of photon-assisted 
Josephson current. We are still curious whether the main resonance located at Eo = 0 can 
also be reversed. In the case of Eo = 0, the current formula is further simplified since | = I 


is required by the resonant tunneling. Fig.3 shows the curves of Io = I(¢ = Z, Eo = 0) vs 


the MW strength o, for four groups of MW frequencies. These numerical results reveal some 
of the properties of Jp: (1) Io has non-monotonous dependence on o, which is natural due 
to the oscillatory nature of Bessel function. (2) For some frequencies, e.g., w — 0.3 in Fig.1, 
lo is always positive; while for others, e.g., w = 0.6, Io reverses its sign in a certain range 
of a. In other words, MW with proper strength and frequency may also induce m-junction 
transition on the main resonance. (3) The dependence of Io on a changes abruptly around 
the frequency of w = A/n (n = 1,2,3---). Roughly, the reason is that when the frequency 


crosses A/n, A+, and B+, jump out of the superconducting gap and feel the singularity in 


the density of states of S electrodes. 

In conclusion, we have shown theoretically that MW irradiation has a non-trivial effect 
on the Josephson current flowing through S-QD-S structure.. The photon-assisted tunneling 
structures appear near Eo = 5hw, and each contains a negative and a positive peak, indi- 
cating m-junction transition induced by MW. The main resonance located at Eo = 0 may 
also reversed with proper MW strength and frequency. The understanding of these results 
requires the knowledge of PAABS, which are formed due to multiple PAAR. We hope that 
this work can stimulate experimental interests in S-QD-S structure, where QD may be a 
nanoparticle or a gate defined geometry in 2DEG. 
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FIGURE CAPTIONS 


Fig. 1 The curves of the critical supercurrent عل‎ vs the resonant energy level Eo, with the 
increase of MW strengths. Parameters are: w = 0.3, T = 0.01. The inset in a = 0 
schematically shows the setup of an S-QD-S structure under the MW irradiation. The 
inset in a = 0.3 demonstrates the formation of PAABS around Eo = —ihw. 


Fig. 2 The analysis of current carrying spectrum j(e) for the point (a)-(e) marked in the 
curves of به‎ = 0.6 in Fig.1. For comparison, the spectrum for رط‎ = 0 (solid) and 


Eo = 0.2 (dotted) without MW irradiation are also shown on the top. 


Fig. 3 The curves of the main resonance Jo vs the MW strength a, for MW frequencies w 


chosen near A/3, A/2, A, and w > A. 
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